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Summary 



We pres ent a practica l way to find the matching priors proposed by 



Welch Sz Peers 



(119631) and 



PeersI (Il965l ). We investigate the use of saddlepoint approximations 
combined with matching priors and obtain p-values of the test of an interest pa- 
rameter in the presence of nuisance parameter. The advantage of our procedure 
is the flexibility of choosing different initial conditions so that one can adjust the 
performance of the test. Two examples have been studied, with coverage verified 
via Monte Carlo simulation. One relates to the ratio of two exponential means, 
and the other relates the logistic regression model. Particularly, we are interested 
in small sample settings. 
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1. Introduction 



We consider inference on a single scalar parameter in the presence of nui- 
sance parameters. Under the frequentist settings, conditional inference can be 
complicated. Bayesian method can simplify frequentist elimination of nuisance 
parameters. The frequentist and the Bayesian approache s can be connec t ed by 



mat ching 



priors . Matching priors were first proposed by 



Welch fc PeersI fll963h 



and lPeersI (1l965l ). Determining a matching prior is equivalent to finding a solution 
of a first order partial differential equation. Only in simple circumstances, such 
as when par ameters are orthog o nal, t he partial differential equation can be solved 



analytically. iLevine fc Casellal (l2003l ) note that "Unfortunately, except for these 
cases, the solution of the resulting partial differential equations becomes quite a 
hurdle; our only hope is to find numerical solutions to these partial differential 
equation." 

We will see a practical way to solve for the m atching priors, without the in- 



volvement of the back transformation described by 



Levine fc Casella (2003). This 



procedure is easy to understand, can be implemented in R, 



R Development Core Team 



( 120071 ) and is suitable to all kinds of initial conditions. 



Our implementation o f matching priors for the approximations proposed by 



DiCiccio fc Martini (1l993l ) is less complicated than other frequentist methods. Di- 
Ciccio and Martin's approximations are saddlepoint approximations that make 
use of Bayesian-frequentist parallels. Our proposed implementation requires less 
computati onal effort compar e d to t he iterative Metropolis-Hasting algorithm de- 



scribed by 



Levine fc Casellal (120031 ). 



We end the introduction with a brief outline of this paper. In §2, We review 
the concepts of matching priors and discuss the circumstance when orthogonal pa- 
rameters are presence. Existing analytical and numerical solutions are reviewed. 
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In §3, we present the procedure for solving matching priors, both analytically 
and numerically Specification of initial condition is discussed. We also provide 
information of R software implementatiori of the solving procedure. In §4, the ap- 



proximations of 



DiCiccio fc Martini (119931 ) are reviewed. The application of using 



matching priors conjuncted with DiCiccio and Martin's approximations is illus- 
trated through examples in section §5. Different initial conditions are specified 
for obtaining various matching priors. Finally, §6 contains the conclusion. 



2. Matching priors 



We consider parametric models with random variables Xi, . . . , X„ having joint 
density function that depends on the unknown parameter vector uo. Suppose uo 
is of length d and uj = {uj^, uj"^, . . . , u'^) = [tp, A) with i/j = uj^, the parameter of 
interest, and the nuisance parameter A = (a;^, . . . , a;'^ ). 



Matching priors were proposed by 



iWelch fc PeersI fll963h and 



PeersI (119651 ). In 



the following, denote the matching prior by vr(-). Let pr^(-|X) be the posterior 
probability measure for ip under prior 7r(-). The upper (1 — a) posterior quantile 
constructed on the basis of a prior density function vr (■?/;) has the property that 
it is also the frequentist limit, such that 

pr < V^(i-")(7r, X)\X} = pr^{^ < ^(^-"^(Tr, X)} = 1 - a + 0{n-^). 



When there are no nuisance parameters, IWelch fc PeersI (1l963l ) showed that 
the appropriate choice of 7r(a;) is7r{uj) oc {i{uj)Y^'^, where i{tu) = E{— d^ /(cj)/ dcj^}, 
and /(■) is the log-likelihood function. In this case, matching priors can be easily 
obtained. 

In the presence of nuisance parameters, iPeersI ( 1l965l ) showed that tt{uj) must 
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be chosen to satisfy the partial differential equation 
d „ d 



d 



4j/-lU-l/2 



} = 0, 



where ijk{oj) = E{—d'^l{u!)/duj^duj'^} and (i-'^) is the dx d inverse matrix of (ijk)- 
If the parameter of interest and the nuisance parameter vector are orthogo- 
nal, solving the partial differential equati on ([H) is rela t ively easy. We follow the 



Cox fc ReidI (119871 ). Orthogonahty is 



definition of parameter orthogonality by 
defined with respect to the expected Fisher information matrix. The most direct 
statistical interpretation of parameter orthogonality is that the relevant compo- 
nents of the original statistic are uncorrelated. In general, it is possible to obtain 
orthogonality of a scalar parameter of interest to a set of nuisance parameters. 

When the parameter of interest ip is orthogonal to a set of nuisance parame- 
ters, equation reduces to 



dip 



dip 



(2) 



Tibshiranil (jl989l ) showed that solutions were of the form 'K{ip^ A) oc {i^^{ip^ •^)}^^^ 
(7(A), where g{X) is arbitrary, and the suggestive notation i^^{ip, A) is used in place 
of iu{ip, A). 

However, choosing a parametrization to achieve parameter orthogonality is 
not always easy, and it can be hard in some cases. It is equivalently hard to ob- 
tain orthogonalization and to solve the partial differential equation ([T]) directly, 
since the orthogonalization procedure also requires solutions to partial differential 
equations of form similar to ([1]). Stai cu and Reid (2007), s tudied the use of match- 



ing priors with the approximation of 



DiCiccio fc Martin 



(I1993I I under orthogonal 



parametrization, and showed that the Peers-Tibshirani class of matching priors 
was essentially unique. One can modify the arguments in this paper to solve the 
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partial differential equation that defines the orthogonality transformation, and 
att empt, using orthogonali ty, to narrow down the class of matching priors. 



Levine fc Casellal (120031 ) proposed a general procedure to solve the partial dif- 
ferential equation ([1]) numerically, in models with a single nuisance parameter. 
Firstly, they transform the parameters into another parameter space, solving 
the equation, and then transform back to the original parameter space. The 
numerical application of this procedure is not necessarily eas y, and usually the 



trans: 



brmation between the two parameter spaces is nontrivial. 



Levine &: Casella 



( I2OO3I ) implemented their procedure using Mathematica. They did not give in- 



structions on initial condition specification, which is a necessa ry component to 



Sweetind ^OQM) 



give specific solution in solving the partial differential equation, 
introduced data-dependent priors that locally approximate the matching priors, 
and his procedure can deal with vector nuisance parameters. 



3. Solving for the matching priors 

In this section, we introduce a procedure to solve the partial differential equa- 
tion in general parametrization. For simplicity, we consider the dimension of 
the parameter space as 2. First, we give analytical form of the solutions, and 
then practical notes will be presented later in this section. 

In the case that d = 2, equation ([1]) is reduced to 

a{ij,\)z^ + b{ij,\)zx = d{ij,\), (3) 

where 

^(^,A) = log{7r(^,A)}, 
aiiJ,X) = {^'\iJ,\)Y/', 
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and 



The coefficient a{ilj,X) is a diagonal element of the inverse matrix of (ijk), so 
a{ip, A) can not be zero. Dividing both sides of by a{ip, A), we have 

^ a{ip, A) a{ip, A) 
This forces the coefficient of to be 1, which simplifies the procedure of finding 
a solution. 

To solve the equation ([T]), it suffices to solve the following ordinary differential 
equations system 

dA b{^,\) dz d{'4j,X) 



d^ 
d7 



(4) 



d s a{ijj, A) ' d s a{ip, A) 
To be more specific with the solution, let us consider the initial conditions 
prescribed along an initial curve /. Suppose that / is given parametrically, in 
terms of a parameter ^, as 

^ = v^(0, A = A(0. 
Then evaluating z{ip, A) at a point on / is equivalent to expressing z as a function 

ofe, 

z = Z{0=z{^{0,A{0}- (5) 

6{^(0,A(0}' 



Here, it is obvious to see that / can not be tangent to the direction 
We then obtain 

^ = ^(s,0, A = A(s,0 
by simultaneously integrating the two equations defined by 

^=1, ^(.0,0 = ^(0, 

as 



a{vl/(0,A(0}. 



dA _ b{^,X) 
ds a('?/'. A) 



(6) 
(7) 



From the third equation in (jl]), the initial condition is given by ([5]). Then we 
have, 

z[sQ,i) = Z[i), (8) 



Equation ([H]) can be integrated by quadrature, once equations (j6]) and ([7]) have 
been solved. 



d{4j{s',0,X{s',0} 



ds'. 



(9) 



These generate a surface in three dimensions, Z{iIj,X), that satisfies both the 
equation ([3]) and the initial condition. When there are no close for m solutions 



for eq uations ([6]), ([7]) and ([8]), numerical solutions can be achieved. 
(1l986l ) presents more mathematical details. 



Rhee et al. 



In obtaining the soluti on formula (|9]) o 



formation as described by 



z(s, we avoid doing back trans- 



Levine &: Casellal (120031 ). Noticing that if we want to 



specify the value of a matching prior at a certain point, say (■?/'*, A*), we can 
directly specify s as ip* and ^ as A* in formula (Q, and then the matching prior 
evaluated at (■?/'*, A*) can be achieved. 

Without loss of generality, set the initial condition 

{vi/(o,A(0,^(0} = (o,e,-i). 

With = 0, we have ip = The equations ([7]) and ([8]) can be simplified as 

dA b{s,X) 



ds a(s. A)' 
dz d{s,X) 
ds a(s. A)' 

We used R package odesolve by 



A(so,0 = A(0, 
z{so,0 = ZiO- 



(10) 



Setzerl (120071 ) to solve equation (ITOl) and get a 



numerical expression of A(-) in s. The command IsodaO in odesolve package is 
designed to solve initial value problems for stiff or non-stiff systems of first order 
ordinary differential equations. It provides an interface to the Fortran ordinary 



different i al equ ation solver of the same name, written by 



Petzold 



HindmarshI (119831 ) and 



(Il983l ). For Qj, we did numerical integration using Sim pson's Rule an d 



CurranI (120051 ) ■ 



employed the R function sintegral () in the Bolstad package by 
Suppose z will be evaluated at (■?/'*, A*). Noticing that A(^) = ^, choose the start 
value as A* in solving (fTOl) . and then choose the upper integration limit as ip* in 
(Q. The procedure is easy to perform if one has an ordinary differential equation 
solver, even if not using the solver provided by R package odesolve. 

Based on the ordinary differential equation iQ, ip = s+\E'(^), i.e. s = tp—'^l^). 
So So must be chosen considering the range of ip. If we choose ^^(0 = 0, then 
ip = s. For the first example in §5, the parameter ip is the ratio of two exponential 
means, and hence ip > 0- Therefore, Sq should be chosen as any positive value. 

In the above we choose the initial values as {"^ , A{^) , Z = (0,,^, — 1). 
Now we will show that the numerical solving procedure is suitable to any initial 
values. 

• Suppose the initial condition for the ordinary differential equation ([7]) is A(so, - 
A(^), for A(^) an arbitrary known function rather than A(^) = ^ as above. The 
solution formula of z is the same as stated in (Q. When solving ([7]), the initial 
value should be chosen as A(A*), no longer A*, if z is evaluated at (V^*, A*). 

• If the initial condition of (EI) is iIj{so,C,) = ^(0; then the solution from the 
equation ([6]) is = s + \E'(^). Therefore, the equation (JTj) becomes, 

dA _ Hg + ^(0,A} 
ds ~ a{s + ^(0, A}' 

Let s = s + \E'(^). By simple change of variables, ([7]) becomes 

dA _ 6(5, A) 
dl ~ 



Equation ([8]) is 



a{s, A) 

dz rf[7A{g-vl/(0,0,A{g-vl/(0,0] 
d~s a[7A{5-vI/(0,O,A{S-vI/(0,O] 
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with z{so — = Z{^), noticing that sq = sq + ^(0- Then the solution of 

z is simply given by the following formula, 

That is to say, the value of the prior on a certain point with the initial condition 
"^("^Oj = ^(0) is obtained by translating the interval of integration when = 
by ^{0- 

• Suppose the initial condition for ([H]) is z(so,0 = ^(0 is a known 

function. This case is even simpler to deal with. One only needs to plug the 
value of Z(0 into (E]). 

Therefore, the suggested numerical solving procedure is suitable to any initial 
values. 

In the above, both the parameter of interest and the nuisance parameter are 
scalars. With dimension 2, it is relatively easy to understand the first order par- 
tial differential equation solving procedure from the geometric point of view, since 
one can draw the initial conditions and the solution surface in a 3-dimensional 



space. In 



Zham 



(120081 ). the solving procedure was extended to multiple nuisance 
parameters, while keeping the parameter of interest as a scalar. The procedure 
of the higher dimension is similar as the one of 2-dimensional model parameters. 
However, when d > 2, it can be computational intensive to implement the proce- 
dure. Also, if there are no explicit expressions for the coefficients in the original 
first order partial differential equation, numerical implementation may be more 
difficult. 



4. DiCiccio AND Martin's approximations 

Likelihood ratio test is widely used in statistical inference. The signed root 
of the likelihood ratio statistic is i? = sgn(-?/' — i/jq)[2{1{lj) — li^ipo, Ao)}]"*^^^, where 
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l{uj) is the log-likelihood function for the unknown parameter vector uj and Aq 
is shorthand for u^^, the constrained maximum likelihood estimator of u. The 
standard normal approximation to the distribution of R typically has error of 
order 0(n~^/^), and R can be used to construct approximate confidence limits 
for ip having coverage error of order 0(n~^^'^). 



Using matching priors, iDiCiccio fc Martin 



(1l993l ) proposed tail probability 



approximations of order 0(n ^). The approximations are sad dlepoint approxima- 



tions that involve Bayesian m ethod. The approx i matio ns of 



DiCiccio &: Martin 



( I1993I ) can be expressed in the iBarndorff-NielsenI (jl980h format 



<^>{R + R-Hog{T/R)}, 



(12) 



(13) 



(14) 



and the iLugannani Sz Ricd (Il980l ) format 

$(i?) + 0(i?)(i?-i-T-i), 

where $ is the standard normal distribution function, and T is defined as 

T = Wo.A„)|-'"W'°-^°'l"'-<r>. 

|-L„WI"''r(V.„,A„) 

Here l^{uj) = dl{u!)/dip, Lc; is the matrix of second-order partial derivatives of 
taken with respect to u; l\x{uj) is the submatrix of l^juji^^) corresponding to 
A; and '/r(c(j) is a matching prior density for u = {ip, A) which satisfies equation 
(Q. Then the resulting approximation is pr(^/' > ipol^) = ^{R + R~^ log(T/i?)}, 
or, pr(?/' > ipo\X) = $(-R) + (f){R){R~^ — T^^). Both of them have relative error 
of order 0(n^^). Approximate confidence limits for ip can be constructed using 
either of (fT2|) or (fT3l) . These confidence limits have coverage errors of order 
0{n^^). To relative error of the order Op{n~^), the variable T is parameterization 
invariant under transforma tions u 1— » ^iP,t(lj)} . 



The approximations of 



DiCiccio fc Martin 



(I1993I ) show their advantages in 



less computational effort compared to the Metropolis-Hastings procedure used 
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by lLevine &: Casellal (120031 ). To calculate T in flMj) . the matching prior requires 
to be evaluated at two points, {ipo, A^q) and a). The initial curve can be chosen 
passing through {ipQ, X^^); that is to say, only the solution on one point u needs 



to be determined. 



5. Examples 
5.1. Ratio of two exponential means 

Let X and Y be exponential random variables with means /i and u respec- 
tively; the ratio of the means, is the parameter of interest. The parameter 
transformation ^/i — Xip~^, v A^/^^j makes the two new parameters il) and A 
orthogonal. Then X and Y have expectations and A?/^^, respectively. 

Suppose we have n independent replications of (X, y). Denote uj = {ip, A). We 
can obtain the log-likelihood function as l{uj) = —n{{ipx + y)/{\y/ilj) + 2 log A} . 

Both approximations of the Barndorff-Nielson format ( |T2l) and the Lugannani 
and Rice format ( |T3l) are considered. Based on these approximations, p-values 
can be calculated. Approximations based on different prior density functions 
mentioned previously may be used to generate an approximate one-sided p- value 
by approximating pr(i? > r), for r the observed value of R. Approximate two- 
sided p- values may be calculated by approximating 2 min{pr(i? > r), pr(i? < r)}. 
One and two-sided hypotheses tests of size a may be constructed by rejecting 
the null hypothesis when the p- value is less than a. Table [1] reports type I error 
probabilities of the 1,000,000 rounds of simulation with n = 10. 

In this example, the parameters ip and A are orthogonal. Using the simplified 
partial differential equation ([21), 7r('?/',A) = l/ip is an explicit solution. Also 
tt{iIj,X) = 1/{iIjX) is another explicit solution. Numerical solutions were also 
calculated. One of the initial condition is {"^ , A{^) , Z = (0,,^, — 1). The 
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resulting matching prior corresponds to the the analytic solution l/ip. Another 
numerically solved matching prior is based on the initial condition (0, S,-, — log.^), 
which corresponds to the the analytic solution !/(■?/' A). From Table [H, one can 
see that the numerical and analytic solutions give almost the same simulation 
results, which confirmed the validity of our numerical solution process. 

Approximations ( fT2l) and (|T3l) have a removable singularity at i? = 0. Conse- 
quently, these and similar formulae require care when evaluating near i? = 0. In 
these cases, for all but the most extreme conditioning events, the resulting con- 
ditional j9-value is large enough as to not imply rejection of the null hypothesis, 
and so these simulated data sets are treated as not implying rejection of the null 
hypothesis. 

5.2. Logistic regression 

We consider a logistic regression model with a binary response Y and only one 
explanatory variable X. Let uoi denote the unknown intercept and uj2 denote the 
unknown effect of the explanatory variable. Suppose uj2 is the parameter of in- 
terest and uji is the nuisance parameter. We will solve matchin g priors and apply 



■ DiCic cio and Martin's approximations to do inference about uj2- 



Levine &: Casella 



( 20031 ) considered a similar example. 

Let Yi be the response variable taking binary values with success proba- 
bility as Pi, and Xj be the explanatory variable following uniform distribution 
?7(0, 1). Suppose there are n independent replications of (Xj,Fj). Fit the model 
\og{pi/ (1 —Pi)} = v[uj = uJi+ u!2Xi, where Vi = (1, Xj)' and UJ2 is the parameter of 
interest. Inverting the equation, we have pj = (1 + e~'"i'^)~^. We can obtain the 
log-likelihood function as l{uj;x) = 1]^=! 2/« logfe/(l ^ Pi)} + Z]r=i log(l - Pi)- 
The first derivative of the log=likelihood function is V'{y —p), where V is the 
design matrix with in row i. The second derivative of log-likelihood function 
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is —V'WV, where is a diagonal matrix with diagonal elements pi{l — pi),i = 
1, ■ ■ ■ ,n. 

Using sample size n = 30, we generate data satisfying the logistic regression 
model with Ui = —1, uj2 = 0.5, and the explanatory variable X following uniform 
distribution [/(0, 1). In this case, generally the parameters ui and UJ2 are not 
orthogonal. We use the numerical procedure described in §3 and study perfor- 
mances of different initial conditions. Table [2] contains type I error probabilities 
for both one-sided and two-sided tests for approximations of both Barndorff- 
Nielson format and Lugannani and Rice format, based on 10,000 rounds of sim- 
ulation. 

As we mentioned previously, approximations (fT2|) and (fT3|) have a removable 
singularity at i? = 0. We deal with this singularity the same way as in §5.1. 

In the following, we give some instructions on how to change the initial con- 
dition and how to choose favorable initial conditions. Initial condition (0,^, — 1) 
gives type I error probabilities larger than the nominal level 0.05; that is to say, it 
has the tendency to underestimate tail probabilities and reject the null hypoth- 
esis. We want to choose initial conditions to obtain a test whose type I error 
rate is closer to the nominal level. We adjust the initial condition when solving 
the partial differential equation ([T]), and use the Barndorff-Nielson format of the 
approximation. The quantity T in (|T4l) is the only part in the approximation 
that relates to matching priors. For a one-sided test, when the probability is 
small and close to 0, R and T are negative. Making -|- log(T/ R)} larger 
is equivalent to making T bigger. Also one may notice that Z{^) is used only 
in equation IQ. Suppose the initial condition is {\E'(^), A(^), Z(^)}. Keep the 
first two components of the initial condition, \E'(0 and A(^), unchanged, and 
only modify the third term, Z[^). By doing so, the integral part in equation 
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(Q is kept unchanged and z varies only with Z{^). By changing Z{C,), we want 
to adjust T to be bigger. Because T is negative when reject a hypothesis, and 
matching priors appear in T as a ratio, one can construct a Z(-) such that the 
ratio, exp{Z(-?/', A)}/ exp{Z{ipQ,\Q)}, will be smaller than 1; recall that 1 is the 
value of the ratio when Z{C,) = —1. Based on the above arguments, Z{-) function 
is constructed as Z(^) = — log{(^ + 1)'^ + 1}, where g is a tuning parameter and 
leads Z{-) to an even function. As an even function, Z{C,) achieves its maximum 
value at —1, where —1 is the true value for the nuisance parameter when data 
were simulated. We have constructed priors using knowledge of the true value 
of the nuisance parameter. Of course, in practice this knowledge is unavailable. 
One might instead use an estimator of the nuisance parameter in place of the 
true value. 

When Z{^) increases quickly, such as g = 2 in table Ej the type I error 
probability deviates far away from the nominal level in the other direction. If a 
more slowly increasing functions is used, the performance of type I error may be 
better. 

Unfortunately, with some choices of initial conditions, such as the last three 
listed in table [21 the Lugannani and Rice format approximation may fall out- 
side the range of and 1 in some cases. For example, the initial condition of 
[0, C,, — log{(.^ + 1)^ + 1}] yielded 5 such probabilities out of 10,000 data sets. We 
convert those values to or 1 by min{max(j9, 0), 1}, where p is the p-value that 
is outside and 1. 

For the parameter of interest 002, we calculate credible intervals using DiCiccio 
and Martin's approximation in Barndorff-Nielson format. With initial condition 
(0, ^, — 1), out of 1,000 generated data sets, there are 938 credible intervals covered 
the true value 0.5. With initial condition [0, 1^, — log{(.^ + 1)^/^ + 1}], for the 
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parameter of interest ijJ2-, there are 954 credible intervals covered the true value 
0.5. 

We apply the above procedure to a real data set from 



Hosmer fc Lemeshow 



(12000 



Table 1.1). The response variable is coronary heart disease indicator, 
and the explanatory variable is age, x. One hundred subjects were included in 
the study; i.e. n = 100. We fit the logistic regression model following the same 
definition as above, with ui defined for the unknown intercept and 002 for the effect 
of age on heart disease status. Using initial condition (0,^,-1) and Barndorff- 
Nielson format approximation, a two-sided testing p-values is 5.532326 x 10~^, 
and five and ninety-five posterior percentiles are of 0.07 and 0.15 respectively. 



6. Conclusion 



Matching priors were first proposed by 



Welch fc PeersI fll963l ) and 



PeersI f )l965l ) 



In the general parametrization, if the parameter of interest and the nuisance 
parameters are not orthogonal, solving the prior from a first order partial differ- 
ential equation is nontrivial. This paper presents a practical way to solve for the 
matching priors and the procedure can be suitable to all ki nds of initial condi- 



tions. Matching priors can be used with the approximations of 



DiCiccio &: Martin 



(I1993I ). By choosing differential initial conditions one is able to improve the per- 



formances of DiCiccio and Martin's approximations. 
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Table 1: Ratio of two exponential means: type I error probability 





BN Format 


LR Format 


Tests 


1-sided 


2-sided 


1-sided 


2-sided 


Likelihood ratio test 


0.0520 


0.0526 


0.0520 


0.0526 


I.e. (o,e,-i) 


0.0456 


0.0441 


0.0456 


0.0441 


Analytic solution: l/-^ 


0.0456 


0.0441 


0.0456 


0.0441 


I.e. (0,e,-logO 


0.0499 


0.0498 


0.0499 


0.0498 


Analytic solution: !/(■?/' A) 


0.0499 


0.0498 


0.0499 


0.0498 



*I.C. stands for initial condition. 

'^Results are based on 1,000,000 rounds of simulation with n = 10. 
*Tests are of nominal type I error 0.05. 



18 



Table 2: Logistic regression: type I error probability 





BN Format 


LR Format 


Test 


1-sided 


2-sided 


1-sided 


2-sided 


Likelihood ratio test 


0.054 


0.060 


0.054 


0.060 


I.e. (o,e,-i) 


0.052 


0.057 


0.052 


0.057 


LC. [O,e,-log{(e+ 1)2 + 1}] 


0.028 


0.019 


0.031 


0.020 


I.e. [o,e,-iog{(e+ 1)2/^ + 1}] 


0.041 


0.041 


0.044 


0.046 


LC. [0,e,-log{(e+l)2/" + l}] 


0.045 


0.048 


0.046 


0.050 



*I.C. stands for initial condition. 

'^Results are based on 10,000 rounds of simulation with n = 30. 
*Tests are of nominal type I error 0.05. 
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